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ABSTRACT 

We present a new Chandra observation of supernova remnant (SNR) G350. 1-0.3. The high resolution X- 
ray data reveal previously unresolved filamentary structures and allow us to perform detailed spectroscopy 
in the diffuse regions of this SNR. Spectral analysis demonstrates that the region of brightest emission is 
dominated by hot, metal-rich ejecta while the ambient material along the perimeter of the ejecta region and 
throughout the remnant's western half is mostly low-temperature, shocked interstellar/circumstellar medium 
(ISM/CSM) with solar-type composition. The data reveal that the emission extends far to the west of the 
ejecta region and imply a lower limit of 6.6 pc on the diameter of the source (at a distance of 4.5 kpc). We 
show that G350. 1-0.3 is likely in the free expansion (ejecta-dominated) stage and calculate an age of 600- 
1200 years. The derived relationship between the shock velocity and the electron/proton temperature ratio is 
found to be entirely consistent with that of other SNRs. We perform spectral fits on the X-ray source XMMU 
J172054. 5-372652, a candidate central compact object (CCO), and find that its spectral properties fall within 
the typical range of other CCOs. We also present archival 24 fim data of G350. 1-0.3 taken with the Spitzer 
Space Telescope during the MIPSGAL galactic survey and find that the infrared and X-ray morphologies are 
well-correlated. These results help to explain this remnant's peculiar asymmetries and shed new light on its 
dynamics and evolution. 

Subject headings: ISM: individual (G350. 1-0.3) — stars; individual (XMMU J 172054.5-372652) — stars; 
neutron — supernova remnants 



1. INTRODUCTION 

Observations of young supernova remnants (SNRs) provide 
our best means of relating the thermodynamic properties of 
shock waves and ejecta distributions to the mechanics of su- 
pernovae and their associated progenitors. Chandra studies 
of young remnants have yielded some of the most spectacular 
advances in SNR research and revealed intriguing and unex- 
pected features of these intricate systems. 

G350. 1-0.3 is a young, luminous SNR in the inner Galaxy. 
Its nonthermal and linearly polarized radio emission led to its 
identification as a SNR by Clark et al. (1973,1975), but a 
4.8 GHz Very Large Array radio image (Salter et al. 1986) 
revealed a bizarre, asymmetric morphology, unlike that of 
the typical shell-like structures. In subsequent SNR catalogs, 
the source was either removed or listed as a SNR candidate 
(Green 1991; Whiteoak & Green 1996). G350.1-0.3 has been 
detected with ROSAT and ASCA, where it was labeled IRXS 
J172106.9-372639 and AX J1721.0-3726, respectively (Vo- 
ges et al. 1999, Sugizaki et al. 2001). An observation with 
XMM-Newton (Gaensler et al. 2008, hereafter G08) revealed 
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the source to be concentrated in a bright clump, coincident 
with the region of strongest radio emission. The X-ray spec- 
trum was found to be well-fit by a two-component model 
(XSPEC model VPSHOCK + RAYMOND), consisting of a 
shocked plasma with electron temperature kTe ss 1 .5 keV, ion- 
ization timescale r « 3 x 10" s cm"-' and large overabun- 
dances of all metals, as well as a collisionally equilibrated, 
low-temperature (kT^ « 0.4 keV) component with solar abun- 
dances. The best-fit value for interstellar absorption was 
Nh « 3.8 X 10" cm-2. G08 proposed that G350. 1-0.3 is in- 
teracting with a molecular cloud seen along its eastern edge 
in the '"CO survey of Bitran et al. (1997). These considera- 
tions as well as data from the Southern Galactic Plane Survey 
(McClure-Griffiths et al. 2005) were used to derive an ap- 
proximate distance of 4.5 kpc. Henceforth we scale quantities 
with d4,5, the distance in units of 4.5 kpc. 

G08 detected an unresolved X-ray source, XMMU 
J172054. 5-372652, to the west of the bright emission and 
proposed it to be a neutron star, likely a central compact object 
(CCO), associated with G350. 1-0.3. Its spectrum was found 
to be a good fit to an absorbed blackbody although a power- 
law model with an unphysically high photon index (F = 5.4) 
provided an acceptable fit as well. G08 carried out a search 
for possible pulsations but found no periodic signals in the 
range of 146 ms to 1 hr 

In this paper, we present a detailed study of G350. 1-0.3 
and its candidate CCO. In § 2 we summarize the Chandra and 
Spitzer observations and data reduction. § 3.1 features a dis- 
cussion of the X-ray and infrared imaging and a comparison 
of various spectral fits. In § 3.2 we use the parameters de- 
rived from the models to analyze the dynamics and evolution 
of G350. 1-0.3. A discussion of the compact source XMMU 
J 172054.5-372652 follows in § 3.3. In § 4 we interpret the 
results of these investigations and attempt to present a coher- 
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Figure 1. Panel a shows a Chandra image of G350. 1—0.3 in the 0.5—8.0 keV 
band. The image has been smoothed by a Gaussian kernel with a cr of 3 pixels 
(1.48 arcsec). The detailed structure of the bright emission region in the east 
is not visible due to the high contrast. The dashed contour surrounding the 
SNR is at a level of 3(t above the background and shows the apparent extent 
of the emission. The white cross indicates the centroid of the contour and is 
regarded as the putative center of the SNR throughout this analysis. Regions 
2—4 were used for spectral extraction (as shown in Figure 3) with parameters 
listed in Table 1. The dashed circle, partly shown in the southwest comer, 
was used to subtract the local background. The aiTow indicates the position 
of the candidate CCO XMMU J172054.5-372652. Panel b shows an enlarged 
image of the bright eastern emission indicated by the square region in panel a. 
Regions la, lb and Ic were used for spectral extraction (as shown in Figure 
2), with parameters for regions la and lb listed in Table 1. The colorbars (in 
units of counts/pixel) indicate biightness levels. 

ent and self-consistent picture of G350. 1-0.3. We conclude 
with § 5, a summary of our findings. 

2. OBSERVATIONS 

In order to obtain a high-resolution map of the SNR emis- 
sion structure, study its spectral properties in unprecedented 
detail and obtain tighter constraints on the spectral fit param- 
eters, we observed G350. 1-0.3 on 2009 May 21 and 22 for 
84 ks with the Advanced CCD Imaging Spectrometer (ACIS- 
S) on board the Chandra X-ray Observatory (observation ID 
10102). Standard data reduction was performed to remove 
hot pixels and flares, resulting in an effective exposure time 
of 83 ks. The 24 ^m and 70 fim images of the spatial 



region corresponding to G350. 1-0.3 were obtained during 
MIPSGAL, an infrared survey of the Galactic plane (Carey 
et al. 2009). The target was observed (in both bands) on 
2006 October 06 with the Multiband Infrared Photometer for 
Spitzer (MIPS) instrument on board the Spitzer Space Tele- 
scope. The raw data were calibrated using the standard BCD 
(basic calibrated data) pipeline to correct for bright pixels 
and flares resulting in exposure times of 2.6 s (24 /im) and 
3.15 s (70 /.tm). The 70 /im data were additionally filtered 
to mitigate instrumental effects (stimflash latents and resid- 
ual background drifts). An 8 /im image was obtained during 
GLIMPSE, an infrared survey of the inner Milky Way Galaxy 
(Churchwell et al. 2009). The target was observed on 2004 
September 07 with the InfraRed Array Camera (IRAC) on 
board the Spitzer Space Telescope. The raw data were cali- 
brated using the standard BCD pipeline to correct for hot pix- 
els and flares resulting in an exposure time of 1.2 s. 

3. ANALYSIS 

3.1. Imaging and Spectroscopy 

A Chandra image of G350. 1-0.3 is shown in Figure la. 
We filtered the data on the 0.5 - 8.0 keV energy band and 
smoothed the image by applying a Gaussian kernel with a a 
of 3 pixels (1.48 arcsec). The source is dominated by an ex- 
tended bright clump of material in the east with fainter emis- 
sion extending far to the west. The dashed contour, shown 
in the figure, surrounds the SNR at a level of 3cr above the 
background and indicates the apparent extent of the emis- 
sion. Throughout our analysis, we adopt a source radius of 
R = 3.3c/4.5 pc (2.5 arcmin) and assume, following the argu- 
ment of G08, that (^4 5 = 1. Although the data reveal an el- 
liptical shell-like structure, a large region in the southwest is 
markedly fainter than the adjacent emission with a flux only 
a factor of ^ 1 .5 above the background. The unresolved X- 
ray source XMMU J172054.5-372652 (indicated by the ar- 
row in Figure la) is seen far to the west of the bright emis- 
sion, significantly displaced from the apparent center of the 
SNR (marked by the white cross in Figure la). Figure lb 
shows the detailed structure of the ejecta region and reveals 
previously unresolved filamentary structures running across 
its center We find that the six regions enumerated in Fig- 
ures la and lb represent the range of spectral characteris- 
tics found in this remnant. The spectra and the correspond- 
ing response files were produced using the specextract script 
in Ciao 4.2. We subtracted the local background spectrum 
from a source-free region of the detector, partly shown as 
the dashed region in the southwest comer of Figure 1 a, and 
grouped the data to a minimum of 15 counts per bin. We 
fit each spectrum by an absorbed nonequilibrium ionization 
(NEI) plane-parallel shock model with variable abundances 
(XSPEC model "VNEI"; Borkowski et al. 2001). Elements 
below Mg and above Fe were fixed at solar abundances since 
their contributions in the fitted bandpass are negligible. The 
uncertainties were calculated at 90% confidence (1.6(t). Al- 
though fixing the absorption at an average global value of 
^ 4.0 X 10^^ cm~^ produced acceptable fits for regions la, lb, 
2, 3 and 4 (the exception is region Ic and is discussed be- 
low), we allow Nh to vary in order to account for possible 
small-scale variations in the absorbing column. The spectra 
and fitted models for regions la, lb and Ic are shown in Fig- 
ure 2 with parameters (for regions la and lb) listed in Table 
1. The fits for regions 2, 3 and 4 are shown in Figure 3, with 
parameters listed in Table 1 . 
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Figure!. Plots of spectra extracted from regions la, lb and Ic overlaid 
with the fitted VNEI models and residuals. The NEI model in XSPEC does 
not include atomic data for the argon emission; a Gaussian component was 
added to the model to account for this feature in the spectral fitting. The fit 
parameters for regions la and lb are hsted in Table 1. 



The spectrum of region la (extracted from one of the 
bright filaments in Figure lb) is roughly representative of 
that found throughout the bright eastern emission. It can 
be fit by an absorbed VNEI model with column density 
Nh ~ 3.8 X 10~~ cm"-, electron temperature fcr^ « 1.4 keV, 
ionization timescale r « 2.2 x 10" s cm~^ and large over- 
abundances of Mg, Si, S, Ca and Fe. These characteristics 
unambiguously demonstrate the presence of hot, metal-rich 
ejecta. The spectrum of region lb (extracted from the other 
bright filament) is similar to that of region la although the 
Mg and Si lines clearly indicate a somewhat lower ionization 
timescale. Although morphologically well-defined, we find 
that the spectral characteristics of the filaments are basically 
consistent with the surrounding regions. 

Although G08 found evidence for variations in elemental 
abundances between nearby regions within the bright east- 
ern emission, our spectral fits of small regions in this vicinity 
show no significant variations within uncertainties. In order to 
compare our results with those of the XMM-Newton observa- 
tion, we attempted to fit the entire ejecta-dominated emission 
(region Ic: the same area as fitted by G08) using a VNEI 
model. As can be seen from the high residuals, this model 
does not provide a good fit, and hence, we do not display the 
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Figure 3. Plots of spectra extracted from regions 2, 3 and 4 overlaid with 
the fitted VNEI models and residuals. The fit parameters are listed in Table 



fit parameters. However, the best-fit parameters are roughly 
consistent with those of regions la and lb, indicating that the 
extremely large number of X-ray photon counts (~ 100,000) 
and small error bars likely magnify small deficiencies in the 
model and result in the poor x' statistics. In addition, it is 
possible that small variations in temperature and elemental 
abundances exist, but are not discernible with the statistics in 
our spectra from small-scale, spatially -resolved regions . Us- 
ing the two-component model VPSHOCK + RAYMOND (the 
same as that used by G08 - their fit results are summarized in 
§ 1) for region Ic also fails to provide a good fit. The incon- 
sistencies between our results and those of G08 are likely ex- 
plained simply by the fact that the long exposure in the Chan- 
dra observation exposes weaknesses in the models that were 
previously hidden due to the fact that the XMM-Newton spec- 
trum has fewer photon counts. 

Region 3 is largely representative of the spectral features 
found along the perimeter of the ejecta region and throughout 
the remnant's western half. The spectrum can be described by 
a VNEI model with absorbing column Nh ~ 4.4 x 10^- cm"^, 
electron temperature kT^ w 0.85 keV, ionization timescale 
T sa 2.7 X 10" s cm"-' and near-solar abundances of all met- 
als. These results suggest an interstellar/circumstellar origin. 
Although otherwise similar, the fit for region 4 indicates a 
somewhat lower temperature and an ionization timescale of 
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Table 1 

VNEI parameters for spectral fits of regions la, lb, 2, 3, 4 



Parameter 


Region la 


Region lb 


Region 2 


Region 3 


Region 4 


Wh" (10^2 cm"^) 


-0.2 


4.0+''3 

-0.1 


3.9!«-3 

-0.2 


4.4^-3 

-0.2 


4.3!«-2 

-0.2 




1 4+0.2 


1 2+''2 


8+<'-' 


9+<' ' 


54+0 04 
"■-'^-0.04 


Mg 


-1.8 


-1.8 


2.6+ii 

-0.7 


-0.6 


-0.2 


Si 




7 9+1.4 
'■^-1.4 




1 6+<'-3 

'■"-0.2 


1 7+o^3 

'■'-0.1 


S 


4 4+4.4 




2.3^;^ 


1 2+<'-2 
'■^-0.2 


5+<'^2 
'-'■^-0.1 


Ca 


6 1+2" 

"■'-1.8 


4 9+2.3 




1 1+'" 
'■'-0.9 


0+4.9 


Fe* 




^•■'-1.4 


(1) 


(1) 


(1) 


t(10" scm-3) 


9 9+0.6 
-0.3 


92+0.10 




9 7+0.5 
'-0.5 


13+0 04 

"■' -0.02 


Flux*^ (10"" ergs cm"^ s"') 


8.9 


8.2 


5.4 


11 


7.2 


Xl/doi 


1.94/157 


2.16/144 


1.89/115 


1.55/178 


2.27/154 



Note. — All uncertainties are statistical errors at 90% confidence, i.6cr. Elements below Mg 

and above Fe were fixed at solai' because their contributions in the fitted bandpass ai'e small 

"The absorption was calculated using the model of Wilms et al. (2000) 

''(1) indicates that the elemental abundance was fixed at solar 

'^The values listed are the unabsorbed fluxes over the energy range 0.5 - 8.0 keV 



r ft! 1 .3 X 10'° s cm"'', lower than that of region 3 by a factor 
of 20. The centroid of the Si Hne in the region 4 spectrum 
is at a slightly lower energy than that of region 3 and since 
the goodness of fit statistics are dominated by the Si line (the 
error bars are smallest in this region), it is this feature that 
causes the drastic difference in the ionization timescale. Re- 
gion 2, adjacent to the bright eastern emission, has somewhat 
enhanced abundances of Mg, Si and S and a highly enhanced 
abundance of Ca. Although similar in temperature to region 
3, the overabundance of metals suggests the possible presence 
of ejecta. 

The model histogram for the spectrum of region 4 slightly 
misses the centroid of the S line complex, suggesting that 
the several bright patches within this region likely vary in 
their ionization states. Fitting the individual clumps, how- 
ever, leads to a similar result with the added handicap of poor 
statistics in the spectra. We thus present the fit to the com- 
posite region with the understanding that a deeper exposure 
is necessary to characterize the spectrum more precisely. The 
spectral fits to regions lb, 2 and 4 do not formally correspond 
to good fits based on the statistics. In the case of region 
lb, this deficiency is due to the fact that the model histogram 
misses the minimum of the Mg line and is too low on the low- 
energy shoulder of the Si line. In the region 2 spectrum, the 
fit is too low on the maximum of the Si line, as can be seen 
from the high residuals in this area. We acknowledge these 
weaknesses in the models and use them to simply make qual- 
itative statements about the approximate distribution of ele- 
mental abundances, temperatures and ionization states, rather 
than using them to derive precise numerical results. 

A 24 fim infrared image of G350. 1-0.3, overlaid with the 
contours from the X-ray data, is shown in Figure 4a. The 
two bright clumps in the east and the south correspond to re- 
gions Ic and 4 respectively. Most notably, the upper filament 
from the X-ray image (region lb) appears to outline precisely 
the eastern edge of the infrared emission while the lower fil- 
ament (region la) is notably absent at 24 /im. An enlarged 
24 um image of the ejecta- dominated region, overlaid with 
the X-ray contours, is shown in Figure 4b. While the angu- 
lar resolution at 24 um is much coarser than that provided 
by Chandra, these images are sufficiently sensitive to iden- 



tify two distinct regions of emission with clear morphological 
similarities to the X-ray emission from the SNR (regions Ic 
and 4 in Figure 1). Figures 4c and 4d show IRAC 8 fim and 
MIPS 70 fim images of the spatial region corresponding to 
G350. 1-0.3. We do not detect any significant emission from 
G350. 1-0.3 at 8 /im and although we see 70 fim emission 
in regions corresponding to the SNR, the emission does not 
differ significantly from regions in the surrounding medium. 
With the current resolution we cannot determine whether or 
not some of the 70 /im emission is associated with the SNR. 
We searched the 3.6 /im, 4.5 jim and 5.8 fim IRAC bands 
but found no apparent emission from G350. 1-0.3 or XMMU 
J172054.5-372652. 

3.2. Dynamics and Evolution 

We analyze the dynamics and evolution of G350. 1-0.3 us- 
ing the numerical study of nonradiative SNRs by Truelove 
and McKee (1999), hereafter TM99. We fit the X-ray spec- 
trum of the SNR's entire western half by an absorbed VNEI 
model in order to derive the average post-shock electron den- 
sity and the mass of the swept-up material. The parameters 
of the model are roughly intermediate to those of regions 3 
and 4. We assume a spherical half-shell with thickness R/12 
for the geometry and a value of 4 for the compression ratio 
of post-shock and pre-shock densities. We also take the ra- 
tio between the electron and atomic hydrogen densities to be 
ne/riH = 1-2, valid for cosmic abundances. With the normal- 
ization parameter from the fit (^ 0.056), we use the relation 

= 0.056 X 10'^ cm-\ (1) 

4-nd- 

where V is the estimated volume and d is the distance, 
to derive an approximate post-shock electron density n^ « 

—1/2 —'X 

5.6d^^ cm and a corresponding swept-up mass of ^ 

-3/2 

2d^ 5 M0, with the obvious caveat that the uncertainties are 
unquantifiable since the visible emission — and hence the 
observed radius — is only a lower constraint on the true 
size of the SNR. G08 derived a substantially higher density 
through two independent methods. However, the first calcu- 
lation assumes ionization equilibrium, inferred from using a 
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Figure 4. Panel a shows a 24 /im image of G350. 1-0.3 taken during the MIPSGAL Galactic survey. The integrated flux of regions Ic and 4 is ~ 2.5 Jy and 
~ 0.5 Jy respectively. Panel b shows an enlarged image of the region corresponding to the ejecta-dominated emission with corresponding X-ray contours at levels 
of 1.1, 2.5, 7.0 and 13 counts/pixel. Panels c and d show 8 urn and 70 /im images of the same spatial region as panel a. The contours in panels a, c and d are 
at levels of 0.55, 1.1, 5.0 and 13 counts/pixels and coiTespond to the X-ray data from Figure la. Regions Ic and 4 correspond to the same regions as in Figure 1. 
The colorbars indicate brightness levels in units of MJy/sr. 



two-component model (which we were unable to fit, as dis- 
cussed in § 3.1). The second calculation assumes a Sedov 
solution — an assumption that is likely not valid for reasons 
discussed below. If the molecular cloud interaction scenario 
(as proposed by G08) is correct, the density throughout the 
ejecta-dominated half of the SNR is likely to be higher than 
that in the western half. Hence, extrapolating our estimate to 

the eastern half of the SNR (for a total of AcT^I'^ M©) is 
likely an underestimate of the true swept-up mass, although 
this figure can be used as a rough lower bound. 

For a progenitor with an ejecta density profile a (typical 
of core-collapse supernovae), TM99 demonstrate that the ra- 
dius and age at which a free-expanding SNR enters the Sedov- 
Taylor (ST) phase are given by 



and 



tsT = 0.132E-]I^Mlfpl"\ 



(2) 



(3) 



where M^j is the ejecta mass and Esn is the explosion energy. 
Here po is the pre-shock mass density and is given by 



1.4 



P0-- 



(1.2)(4) 



mnrie = 0.30m//«(. 



(4) 



where nin is the mass of the hydrogen atom, 1 .4 is the equiv- 
alent molecular weight of the hydrogen and heUum mix- 



ture, (assuming cosmic abundances) and the factors 1.2 and 
4 are the electron/hydrogen density ratio and post-shock/pre- 
shock compression ratio, respectively. Using typical values 
of Mej =4-14 Mq for the ejecta mass consistent with the 
formation of a neutron star (Woosley et al. 2002), we calcu- 
late RsT to be 4.1-6.2 pc, somewhat larger than the observed 
radius for G350. 1-0.3 {R = 3.3ii4.5 pc). Assuming that the ob- 
served radius of the SNR is approximately equal to its true 
radius (in § 4 we discuss the scenario in which this assump- 
tion is not valid), this result suggests that G350. 1-0.3 is likely 
in the free expansion (ejecta-dominated) stage and has not yet 
entered the ST phase. 

If we make an assumption about the explosion energy, we 
can determine the age using the relation 



t = 0.988fsr ( 

RsT 



7/4 



(5) 



as derived from TM99. We explore the parameter space for 
Mej = 4-14M0 and £5/^ = 0.5- 1.0 X 10^' ergs and derive an 
age of f = 600- 1200 years. The shock velocity and proton 
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Table 2 

Derived Parameters for the Evolution of G350. 1—0.3 



Radius (arcsec) 

Figure 5. Azimuthally-averaged radial profile of CCO XMMU 
J172054. 5-372652 in the 0.5 — 8.0 keV energy range. The histogram with 
the vertical error bars represents the data while the curve i s the model point 
spread function (PSF). 



temperature follow frortiQ 



and 



kTp = 0.nmHvl 



t 

tsT 



-3/7 



(6) 



(7) 



The results of these calculations are summarized in Table 
2. Using kTe w 0.7 keV for the electron temperature of the 
ISM component, the average value from our fits, the inferred 
electron-to-proton temperature ratio is 0.07-0.27. This is 
in excellent agreement with measured values of this ratio for 
SNR shocks in the inferred velocity range (Ghavamian et al. 
2007). 

3.3. The Compact Source 

Using the derived age for the remnant, we calculate the pro- 
jected velocity of XMMU J172054.5-372652 to be 1400- 
2600 km s~', assuming a displacement of 1.7^^4.5 pc (1.3 ar- 
cmin) from the adopted center of the SNR. We generated a 
lightcurve of XMMU J172054.5-372652 but found no short- 
term flux variation within the Chandra exposure. We ap- 
plied barycenter corrections to the photon arrival times and 
searched for periodicity using the Z^-test (Buccheri et al. 
1983). No pulsations with period longer than 6.4 s and pulsed 
fraction larger than 16% were detected at 99% confidence. 
The Chandra image in Figure la shows no extended emission 
surrounding the CCO and a generated radial profile (Figure 5) 
of the source confirms that it is fully consistent with a model 
PSF. As a note, photon pile-up is negligible with a level of 
less than 5%. 

* This expression for the shock velocity comes from TM99 Table 7. It 
should be noted, however, that the authors overlooked the minus sign in the 
exponent. 



Parameter 




Value 


Age, t, (years) 


600 


-1200 


Shock Velocity, I'j, (kms"') 


1500 


-2900 


Proton Temperature, kTp, (keV) 




3-10 


Speed of CCO, v„o, (kms"') 


1400 


-2600 



We extracted the ~ 4000 CCO counts from a 2 arc sec- 
radius aperture and grouped the source spectrum to a mini- 
mum of 25 counts per bin. Using the SHERPA environment, 
we fit the spectrum to simple absorbed blackbody (BB) and 
power-law (PL) models. Although the PL fit is slightly bet- 
ter than the BB, the large photon index (F = 5.5) suggests a 
thermal origin for the emission. Both fits are consistent with 
those discussed in G08. Since CCOs are often characterized 
by two-component models (BB-i-PL and BB-i-BB), we inves- 
tigated those as well. The results are summarized in Table 
3. The best-fit absorption column density is ~ 4 x 10^^ cm- 
and agrees with that of the SNR. While the CCO spectrum 
is dominated by a BB with kT « 0.4 keV and a small emis- 
sion radius of ~ 3 km, the high energy component is not very 
well constrained. Our results suggest that it can be modeled 
by either an additional BB component with a very high tem- 
perature (0.9 keV), or a PL component. Adding a second 
component improves the fit although not at a statistically sig- 
nificant level. Unfortunately, the current data do not allow 
us to distinguish between these two scenarios and a deeper 
exposure is required. We present the spectrum of XMMU 
J 172054.5-372652 overlaid with the BBh-PL model in Figure 
6. 

4. DISCUSSION 

The high resolution Chandra data reveal that the emission 
from G350. 1-0.3 extends far to the west of the bright ejecta 
region and provide a rough constraint on the size of the SNR, 
with the caveat that the asymmetric morphology makes it dif- 
ficult to extract a reliable diameter G08 estimated the age of 
G350. 1-0.3 to be ^ 900 years based on the diameter of the 
ejecta region alone {D w 2.6dm pc) and the shock velocity 
derived from the electron temperature. However, we see that 
the emission from G350. 1-0.3 extends far beyond the bright 
ejecta-dominated region and our analysis (using the model 
treatment of TM99) shows that the electron temperature is 
almost certainly an underestimate of the proton temperature. 
With the diameter, a post-shock density estimate and an as- 
sumed range of values for M^j and Esn as input parameters, 
we calculate a shock velocity of 1500-2900 km s"' , a proton 
temperature of 3 - 10 keV and an age of 600- 1200 years. We 
see that although the size of the SNR is far larger, the age is 
consistent with earlier estimates due to the lower ISM density 
estimate. 

The enhanced elemental abundances and high temperature 
of the bright emission in the remnant's eastern half unambigu- 
ously demonstrate the presence of hot, metal-rich ejecta. By 
contrast, spectral analysis of the emission around the perime- 
ter of the ejecta region and throughout the SNR's western 
half (with the possible exception of region 2 — see discus- 
sion in § 3.1) reveals cooler material with solar-type com- 
position, which supports shocked interstellar/circumstellar 
medium (ISM/CSM) for its origin. This well-defined inter- 



Chandra Observation of G350. 1-0.3 7 



Table 3 

Spectral Fits to XMMU J172054.5-372652 



Pai'ameter 


BB 


PL 


BB+PL 


BB+BB 






4.4 ±0.1 






kTi (keV) 


0.50 ± 0.01 




A /n+0 05 

o.42;[;;^5 


0.41 ±0.05 


R, (km) 


1.59 ±0.02 




2.59 ±0.06 


2.94 ±0.03 


r 




5.5±0.3 


3+- 




kT2 (keV) 








'-'•'-0.2 


R2 (km) 








0.19±0.01 


Flux" (10"'^ ergs cm"^ j-l) 


4.4±0.1 


4.7 ±0.2 


4.9 ±0.2 


4.7 ±0.1 


Fi/F2^ 






1.5 


9.8 


xl/doi 


1.02/105 


0.94/105 


0.86/103 


0.86/103 



"The flux is over the energy range 0.5 - 10.0 keV and has been coixected for fore- 
ground absorption 

''Unabsorbed flux ratio between the first and second spectral components in the 1.0— 
10.0 keV energy range 



10 



(BB+PL) 

1 













energy (keV) 

Figure 6. The X-ray spectrum of CCO XMMU J172054.5-372652 overlaid 
with a BB-hPL model. The fit parameters are shown in Table 3. 

face between the ejecta and shocked ISM regions further sup- 
ports G08's suggestion that G350. 1-0.3 is interacting with a 
molecular cloud in the east. If this is indeed the case, the ex- 
panding shock wave propagating eastward through the dense 
gas would form a strong reverse shock and heat the metal-rich 
ejecta to X-ray emitting temperatures while the shock front 
propagating in the other direction would still be expanding 
relatively unimpeded into the diffuse ISM. 

Although we cannot conclusively determine whether the 
24 i^im emission is due to dust or line emission, the simi- 
larity between the infrared and X-ray morphologies supports 
the shocked dust scenario (Seok et al. 2008). In addition, the 
lack of significant emission at 70 fim and in the IRAC band is 
consistent with dust emission that peaks between those bands. 
The catalog of 24 fim MIPSGAL sources by Mizuno et al. 
(2010) lists only a handful of SNRs, but spectroscopic obser- 
vations with Spitzer have shown that SNRs interacting with 
dense molecular clouds produce a number of shocked molec- 
ular and atomic species in the mid-infrared band (Hewitt et 
al. 2009). Thus the presence of 24 /im emission coinciding 
with the eastern bright X-ray clump in G350. 1-0.3 is con- 
sistent with a molecular cloud interaction. As can be seen 
from Figure 4, the upper filament (region lb) appears to out- 
line precisely the eastern edge of of the region corresponding 
to the ejecta-dominated emission. By contrast, the lower fil- 
ament (region la) is not visible at 24 ^m. The Mg and Si 
lines in the spectra of these two regions (Figure 2) as well 
as the spectral fit results (Table 1) suggest that the ionization 



timescale of the region la filament is noticeably higher than 
that of the region lb filament. Although we cannot determine 
whether the difference is due to a higher density or longer 
time since shock, the material in region la (in either case) 
may have cooled/sputtered more than that in the region lb fila- 
ment. Incidentally, both regions of emission at 24 /im (Ic and 
4) coincide with the regions of lowest ionization timescale, as 
determined from their X-ray spectra. This result seems con- 
sistent with the scenario that the 24 jiva emission is due to 
dust heated by the reverse shock propagating back into the 
SNR and the difference in brightness between the various re- 
gions could be attributed to dust cooling or sputtering. 

If the explosion occurred approximately at the center of 
the X-ray emission (RA. 17:21:01.019, DEC. -37:26:49.61), 
the derived age puts the projected velocity of XMMU 
J172054.5-372652 at 1400- 2600 km s"', assuming it was 
formed in the same explosion. The lower limit of this veloc- 
ity, although high, falls within the 1100- 1600 km s~' mea- 
sured for the CCO in Puppis A (Hui & Becker 2006: Win- 
kler & Petre 2007). The upper limit is significantly higher 
than that measured for any CCO. Although the spectrum of 
XMMU J172054.5-372652 can be fitted with a BBh-BB or a 
BB-i-PL model, the former requires a BB component of unusu- 
ally high temperature. On the other hand, a BB-i-PL model is 
not uncommon when comparing to the X-ray spectra of other 
CCOs. The best-fit BB-i-PL parameters in Table 3 fall within 
the typical range of others (see Table 2 in Pavlov et al. 2004). 
Indeed, XMMU J172054.5-372652 could be a close cousin 
of the CCOs inside SNRs Vela Jr (G266. 1-1.2) and Puppis A. 
These objects have an age of 1 - 3 kyr and a blackbody tem- 
perature of 0.4 keV. We also note that for the CCO inside 
Puppis A, XMM-Newton observations indicate an additional 
hard spectral component that can be fitted with either a PL of 
r = 2.0-2.7 or a hard BB with kT = 0.5-1.1 keV (Becker & 
Aschanbach 2002), very similar to our case. 

In addition to the prominent difference in flux between the 
bright eastern emission and the surrounding regions, there 
is an intriguing contrast in brightness between the ambient 
material in the far west and the bright regions of shocked 
CSM/ISM (regions 3 and 4 in Figure la). We estimate the 
brightness ratio of these regions to be ^ 3, implying a density 
ratio of ~ 1 .7. Thus, if more of the remnant lies unseen fur- 
ther to the west, it would need to be only slightly less dense 
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to emit below the background. If this is indeed the case, it 
would change all the derived SNR parameters and imply that 
the SNR may have already entered the ST phase. The larger 
radius would increase the SNR's age, lower the shock velocity 
and lower the proton temperature. In addition, it would im- 
ply that XMMU J172054.5-372652 may be closer to the ex- 
plosion center than evident from the visible emission, which 
would place its speed at a lower, more typical value. On the 
other hand, if G350. 1-0.3 is indeed interacting with a molec- 
ular cloud along its eastern edge, the apparent center of the 
emission is likely a poor estimate of the explosion location. 
In this case, the supernova center would be closer to the east- 
em edge, which would tend to increase the velocity of the 

ceo. 



5. CONCLUSION 

In this paper, we have investigated the dynamical prop- 
erties of G350. 1-0.3 and its candidate neutron stai" XMMU 
J172054. 5-372652. By doing spectral modeling on the high 
resolution X-ray data, we have found that the region of bright- 
est emission is dominated by hot, metal-rich ejecta while the 
diffuse material throughout the surrounding regions is mostly 
cooler swept-up CSM/ISM with solar abundances. These re- 
sults further support the conclusion of Gaensler et al. (2008) 
that G350. 1-0.3 is interacting with a dense molecular cloud 
in the east. The X-ray imaging has resolved new morpho- 
logical features and revealed that the SNR is far more ex- 
tended than apparent from the ejecta region alone. We used 
the numerical model of Truelove and McKee (1999) to show 
that G350. 1-0.3 is likely in the free expansion phase and de- 
rive values for its age, shock velocity and proton temperature. 
We examined the relationship between the inferred shock ve- 
locity and the derived electron/proton temperature ratio for 
G350. 1-0.3 and found it to be entirely consistent with that 
of other SNRs. The derived age puts the speed of XMMU 
J172054. 5-372652 at an unusually high value although its 
spectral characteristics are found to be consistent with CCOs 
in other remnants. We presented an archival 24 /im image 
of G350. 1-0.3 and found that the infrared and X-ray mor- 
phologies are well-correlated, offering additional support for 



a molecular cloud interaction. 



The authors would like to thank Tea Temim and Daniel Cas- 
tro for their input and several helpful discussions. I.L. ac- 
knowledges support from Chandra grant GO9-0059X. RO.S. 
acknowledges partial support from NASA contract NAS8- 
03060. J.D.G is supported by an NSF Astronomy and Astro- 
physics Postdoctoral Fellowship under award AST-0702957. 
B.M.G. acknowledges the support of the Australian Research 
Council through grant FF0561298. 

REFERENCES 



Becker, W., & Aschenbach, B. 2002, in Proc. 270, WE-Heraeus Seminar on 
Neutron Stars, Pulsars, and Supernova Remnants, ed. W. Becker, H. 
Lescli, & J. Tramper (MPE Report 278; Garching: MPI), 64 

Bitran, M., Alvarez, H., Bronfman, L., May, J., & Thaddeus, P 1997, 
A&AS, 125, 99 

Borkowski. K. J., Lyerly, W. J., & Reynolds, S. P 2001, ApJ, 548, 820 
Buccheri, R., et al. 1983, A&A, 128, 245 
Carey, S. J., et al. 2009, PASP 121, 76 
Churchwell, E., et al. 2009, PASP 121, 213 
Clark, D. H., Caswell, J. L., & Green, A. J. 1973, Nature, 246, 28 
Clark, D. H.. Caswell, J. L., & Green, A. J. 1975, Aust. J. Phys. Astrophys. 
Suppl., 1 

Gaensler, B. M., et al. 2008, ApJ, 680, L37 (G08) 

Ghavamian, P, Laming, J. M.. & Rakowski, C. E., 2007, ApJ, 654, L69 

Green, D. A. 1991, PASP 103, 209 

Hewitt, J. W., Rho, J., Andersen, M., & Reach, W. T. 2009, ApJ, 694, 1266 
Hui, C. Y., & Becker, W. 2006, A&A, 457, L33 
McClure-Griffiths, N. M., et al. 2005, ApJS, 158, 178 
Mizuno, D. R., et al. 2010, AJ, 139, 1542 

Pavlov, G. G., Sanwal, D.. & Teter, M. A. 2004, in lAU Symp. 218, Young 
Neutron Stars and Their Environments, ed. F. Camilo & B. M. Gaensler 
(San Francisco, CA: ASP), 239 

Salter, C. J., Patnaik, A. R., Shaver, P A., & Hunt, G. C. 1986. A&A, 162, 
217 

Seok, J.Y, Koo, B.-C, Onaka, T., et al. 2008, PASJ, 60, 453 

Sugizaki, M., Mitsuda, K., Kaneda, H., Matsuzaki, K., Yamauchi, S., & 

Koyama, K. 2001, ApJS, 134, 77 
Truelove, J. K., & McKee, C. F 1999, ApJS, 120, 299 (TM99) 
Voges, W., et al. 1999, A&A, 349, 389 
Whiteoak, J. B. Z., & Green, A. J. 1996, A&AS, 118, 329 
Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 
Winkler, R F, & Petre, R. 2007, ApJ, 670, 635 

Woosley S. E., Heger A., & Weaver T. A., 2002, Rev. Modem Phys., 74, 
1015 



